Electronic Structure near Impurities in the Chains of YBa 2 Cu 3 06+i 
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We study the electronic structure near impurities in the one- dimensional chains of YBa2Cu306+a;. 
Assuming chain superconductivity below T c due to a proximity coupling, we show that only a 
magnetic impurity induces resonance states in the local density of states. Its spatial structure 
reflects the particle-hole nature of chain superconductivity and thus distinguishes it from other 
broken symmetry phases. For a sufficiently large scattering strength, or due to interference effects 
arising from a second impurity, the chains undergo a quantum phase transition to a polarized state. 

PACS numbers: 72.10.Fk, 71.55.-i, 74.50. +r, 
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Studying the electronic structure of the CuC>2 planes 
in YBa2Cu306+a; (YBCO), an important step towards 
our understanding of superconductivity, has been com- 
plicated by the presence of one dimensional (ID) CuO- 
chains in this material. The coupling of the CuO-chains 
to the CuC>2-planes changes the electronic structure of 
the latter and might facilitate the occurrence of other 
ID structures, e.g., stripes, in the planes. Much theo- 
retical work has focused on this coupling and its 
implications for a series of experimental probes ranging 
from inelastic neutron scattering || , nuclear quadrupole 
resonance (NQR) infrared spectroscopy Q, angle- 

resolved photoemission |tJ and scanning tunneling mi- 
croscopy (STM) experiments ||. In turn, the same cou- 
pling necessarily changes the chain electronic structure 
when the CuC>2 planes become superconducting (SC), as 
was recently evidenced by NQR and STM Q exper- 
iments. In particular, NQR experiments observe an in- 
crease in the 63 Cu NQR line width near the SC transition 
H , while STM experiments |J find evidence for large os- 
cillations in the local density-of-states (DOS). Whether 
these observations arise from the presence of a charge- 
density wave (CDW) 0), or from impurities in SC chains 
is currently a topic of intense discussion. 

The study of the electronic structure near impurities in 
2D superconductors [ p~0|JTl] ] has been a successful method 
to distinguish between various broken symmetry ground 
states. In this Letter we extend these works to study 
the effects of chain impurities, both magnetic and non- 
magnetic, on the local DOS assuming that the chains 
become SC due to a proximity coupling. Both types of 
impurities induce large spatial oscillations in the DOS, 
and thus the electron density, whose amplitude increases 
below T c . However, in contrast to a non-magnetic impu- 
rity, a magnetic impurity induces two resonance states in 
the chain DOS. Moreover, we present a prediction for the 
spatial dependence of these resonances which reflects the 
particle-hole nature of chain superconductivity. This pre- 
diction therefore discriminates between oscillations aris- 
ing from a CDW and those due to impurities. Finally, we 
find that as the distance between two magnetic impuri- 
ties is reduced, the resonance states split and the system 
can be driven through a quantum phase transition. Our 
work is complementary to that of Eggert |L2] who studied 



the effects of impurities in a ID Luttinger liquid. 

Since in the following we employ a model described in 
detail in Ref. jlll], we only briefly review it here. The 
mean-field BCS-Hamiltonian of the chain-plane system 
is given by |Ul|] 
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where cl, d k are the fermionic creation operators in the 
plane and chain, respectively. One has for the planar and 
chain tight- binding dispersions, et and (k, 

ek = — 2i p ^cos(fc x ) + cos(fc v )J — 4:t' p cos(k x ) cos(k y ) — \x v , 
Ck = -2i c cos(A; x ) - [i c . (2) 



Here, t p = 300 meV, t' p 



0.812t p are the 
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planar hopping elements and the chemical potential [Q, 
and Ak = A (cos(k x ) — cos(fc y ))/2 is the planar SC gap 
with Ao ~ 40 meV |14| ]. The chain Fermi momentum, 
k% ps 0.257T, HH yields fi c = — 1.4Lt c , and we take t± = 
OAtp p3| . After integrating out the planar fermions, the 

chain Greens function G = -(^*/w( T )*iLro( T ')) witn 



= (4.,i,T'4.,i,i> d -*-.i.T» d -fc.M) 



and m, I chain indices, is given by 

q-1 _ ( B(k x ,u> n )a C{k Xl uj n )aa 
\C(k x ,uj n )a -B(-k x ,-u) n )a 

where aa is the unit matrix in spin space, and 

B = (iuj n - Ck) 8l,m 

iky (l — m) 
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In contrast to previous work Q, we showed |H| that 
the experimentally measured clean DOS can only be ob- 
tained qualitatively as well as quantitatively by making 
the ansatz 
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FIG. 1. Chain DOS in the SC state for the clean case 
(dashed line) and at the site of a magnetic impurity with 
/3 = 200 meV and U a = (solid line). We set the lattice 
constant an = 1. Inset: Fermi sea of the chain (solid line) and 
Fermi surface of the plane (dotted line). 
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FIG. 2. Chain DOS as a function of distance from the im- 
purity site for the resonance states at ui res — ±2.6 meV and 
(see upper inset) for the coherence peaks at the gap edges, 
^>c P = ±16.0 meV. A gray-scale intensity plot of the DOS at 
±ujres is shown in the lower right. 



G c 1 (k x ,l,m,Lu n ) = G c 1 (k x ,Lu n ) 5 tt% 
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We argued that the physical origin of Eq.([|), which de- 
scribes the absence of coherent correlations between the 
chains, or equivalently, the coupling of a single chain, lies 
in the presence of planar stripes and/or doping inhomo- 
geneities in the chains. 

We employ the T-matrix formalism |l5) to obtain the 
local DOS 

N(u J ) = --lm\G^ 1 \r,r,uj + ) + G ( - 22 \r,r 7 u J+ )\ , (7) 

7T L J 

in the vicinity of an impurity, with oj + = u> + i0 + . Here, 
one has for the full chain Greens function 



G(r, r', iw n ) = G (r- r', iu n ) + 

G (r, ioj n )f(iu> n )G (-r', iu n ) 



(8) 



where T(iu) n ) represents all scattering processes and Go 
is obtained from Eqs.(^) and (^). In what follows we con- 
sider quite generally an impurity which is both a mag- 
netic and potential scatterer. Assuming for simplicity 
that the spin of the impurity is aligned along the z-axis 
pa we obtain for the scattering Hamiltonian 



(9) 



with V — T3 (Uo<Jo — J Sz&s) and Ti,ai are the Pauli- 
matrices in Nambu and spin-space, respectively. Uq (J) 
is the potential (magnetic) scattering strength and S z is 
the z-projection of the impurity spin. One then finds 



T(itu n ) = 1 - VGo(r x = 0, iuj n ) 



V 



(10) 



We remind that an impurity bound state exists at a fre- 
quency ujres, if det[r _1 (w res )] = 0, while for a resonance 
state only the real part of det[T _1 ] vanishes. 

We first study the case of a purely magnetic impurity 
with Uq = and J / 0. A straightforward analysis 
of Eq.([l(]) shows that there exist two resonance states 
whose frequencies, ±oj res , are given by the solution of 
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where 7 = Vp//3, Vp is the chain Fermi velocity and (3 — 
JS z /2. Moreover, a = t\/v, v = (de^/dky) Po , where 
Po is the momentum on the planar Fermi surface (FS) 
with pg = k F (see inset of Fig. [l]) and A = A(p ). While 
in the limit 7 <C 1, the resonance frequencies are located 
close to ±A, uj res shifts to lower energies with increasing 
7. At the critical 7 C = 1 (i.e, j3 c — vp) the resonance 
states cross at u: res — and the magnetic impurity forms 
a bound state with a spin-up (down) quasi-particle for 
J < (> 0). Accordingly, the ground state of the system 
undergoes a quantum phase transition and changes from 
an unpolarized one for 7 < j c to a spin-polarized one with 
total spin S — 1/2 for 7 > -f c . A similar transition also 
occurs in superconductors with s-wave symmetry |l7| , |l8| . 

In Fig. [l] we present the SC chain DOS, obtained from 
a numerical evaluation of Eq.(||), in the clean case and at 
the site of a magnetic impurity with — 200 meV and 
Uq = 0. A comparison with the clean DOS shows the 
emergence of two resonance states at the impurity site 
with uj res w ±2.6 meV, in agreement with the discussion 
above. Simultaneously, the impurity suppresses the co- 
herence peaks at w cp = ±16.0 meV which are present in 
the clean DOS. Moreover, the impurity induces spatial 
oscillations in the chain DOS, which we present in Fig. 
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for the two resonance states and the coherence peaks at 
the gap edges (see inset). The oscillations at ±uj res pos- 
sess a broad range of wave- vectors with < k osc < 1.45, 
which gives rise to the rapid decay of the oscillations' 
amplitude with distance from the impurity. The decay 
length, £ « 10 — 15ao is consistent with the experimen- 
tally obtained value of AO A ||. Note that the oscilla- 
tions at ±w res are out-of-phase, such that a maximum in 
the DOS for +Lu res coincides with a minimum for — Lu res , 
as is particularly apparent when one plots the spectral 
weight of the resonance states at their peak positions 
(see inset (b)). The complementary intensity pattern 
in real space of the hole-like and particle-like resonance 
states is a reflection of the particle-hole nature of Bo- 
goliubov quasi-particles jll|,^9|. This characteristic fea- 
ture of chain superconductivity is qualitatively different 
from the response of a CDW, and hence discriminates be- 
tween proposed broken symmetry phases in the chains. 
In contrast, the oscillations at ±w cp , the frequency of 
the coherence peaks (see inset), are in-phase. Moreover, 
their amplitude decays much slower with distance from 
the impurity site since the wave-vectors of these oscilla- 
tions are centered in a narrow range, Ak 0.02, around 
k sc = 0.81 (see arrow in the inset of Fig. |l]). The wave- 
length, A, of these oscillations exhibits an asymmetric 
frequency dependence with A(— loq) > X(loq) for u>q > 
and a minimum at uj « 0, in agreement with STM exper- 
iments J^]. The spatial dependence of the DOS oscilla- 
tions is qualitatively robust against changes in the form 
of the planar FS. However, their amplitude increases in 
the SC state, and accordingly, the spatial electron density 
changes. While this result is consistent with the broad- 
ening of the NQR line in the SC state , a quantitative 
analysis requires a microscopic model for the coupling of 
the electron density to the quadrupolar moment which is 
beyond the scope of this Letter. 

We next consider a purely non-magnetic impurity with 
J = and Uq ^ 0. We remind that a potential scatterer 
only induces a resonance state if the SC phase varies 
along the FS, or between Fermi points. However, the 
two fc^-points in the chain possess the same SC phase, 
which leads to dct[r _1 ] ~ 1 + (Uo/v F ) 2 , and no reso- 
nance states exist, in analogy to the s-wave superconduc- 
tor. Our numerical evaluation supports this conclusion 
revealing no resonances but again large spatial oscilla- 
tions in the DOS. 

Since an impurity is in general a magnetic and poten- 
tial scatterer we next study the case J, Uo ^ and present 
in Fig. U the chain DOS at the impurity site for (3 = 200 
meV and three values of Uo. We find that with increas- 
ing Uo, the position of the resonances shifts up in energy 
from tcVes = 2.6 meV for Uo = to u> res = 9.5 meV for 
Uo = 200 meV. This implies, that the critical value j c 
(and thus (3 C ) where the ground state of the system be- 
comes spin-polarized increases , ag ain similar to the case 
of an s-wave superconductor |18fl . Moreover, while the 
spatially integrated spectral weight of the particle-like 
resonance is 1.43 times larger than that of the hole- like 
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FIG. 3. Frequency dependence of the chain DOS at the im- 
purity site for (3 = 200 meV and U = (solid line), U = 100 
meV (dashed line), and Uo — 200 meV (dotted line). 



resonance for Uo — 0, the ratio is 0.77 for Uo = 200 meV. 
Thus spectral weight is transferred from the particle-like 
to the hole-like resonance with increasing Uo, however, 
their phase-shifts remain unchanged. 

Due to the ID character of the chains, interference ef- 
fects between two impurities should be large and signifi- 
cantly change the single-impurity DOS discussed above. 
To study these effects, we employ a T-matrix formalism 
for two identical impurities where 

G(r, r', iu> n ) = G (r, r', iui n ) + 

2 

^ G (r, T m ,iu n )f(T iw n )G (r„,r' ,iui n ) , (12) 

n,m— 1 



ri : 2 are the positions of the impurities, and T(r„ 
obeys the Bethe-Salpeter equation 

T(r m ,r ni iu n ) = VS rm .r n 



+Vy^G (r m , n, iu) n )f(Y h Y n ,iun) ■ (13) 
1=1 

In general, we expect that the chain DOS for two impu- 
rities is similar to that of a single impurity if the inter- 
impurity distance, A7' = ri— r±, is larger than the spatial 
extent of single impurity effects, £ « 10 — 15 at Lu res (see 
Fig. ||). In contrast, for Ar < £ significant interference 
effects from multiple scattering processes should occur. 
To test this conjecture, we present in Fig. ^ the numeri- 
cally obtained chain DOS at one of the two impurity sites 
as a function of Ar. For Ar = 00 (Fig. [|a, solid line) we 
recover as expected the DOS for a single impurity with 
two resonance states at uj res = ±2.6 meV (see Fig. [l]). 
As the distance between the impurities decreases, the 
interference effects arising from multiple scattering pro- 
cesses lead to the formation of bonding and anti-bonding 
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FIG. 4. Frequency dependence of the chain DOS for two 
identical impurities at one of the impurity sites with f3 = 200 
meV, Uo = and (a) Ar = oo (solid line), Ar = 20 (dashed 
line), and Ar = 15 (dotted line), and (b) Ar = 10 (solid line), 
Ar = 5 (dashed line), and Ar = (dotted line). 



resonances state, and thus to a splitting of the two res- 
onances into four (Ar = 20, dashed line). This result 
provides a possible explanation for the experimental ob- 
servation of four resonance states in the chain DOS ||. 
The splitting first becomes discernible when Ar « 2£, in 
agreement with our conjecture above. The splitting in- 
creases until the inter-impurity distance reaches Ar = 15 
(dotted line). For even smaller Ar = 10 (Fig. ||b, solid 
line), the peaks recombine again and spectral weight is 
transferred to new resonance states at higher energies 
(Ar = 5, dashed line). These new states are precursors 
of the resonances one obtains in the limit Ar — > (dot- 
ted line) when the DOS is that of a single impurity with 
twice the scattering strength, j3 = 400 meV. Note that 
(3 = 200 meV is below the critical value, (3 C , to induce 
a transition of the ground state, while (3 — 400 meV is 
above this critical value. Thus, the system can be driven 
through a quantum phase transition from an unpolarized 
to a polarized state by continually varying the distance 
between two or more impurities. 

In conclusion we have studied the electronic structure 
near impurities in the ID chains of YB^Cvl^Oq+x. As- 
suming that the chains are SC below T c due to a prox- 
imity coupling, we find that a magnetic impurity induces 
two resonance states in the local DOS. The spatial oscil- 
lations of these resonances reflect the particle-hole nature 
of the chain superconductivity, and are therefore charac- 
teristic features which distinguish it from other broken 
symmetry phases. The amplitude of these oscillations 
increases below T c , which is consistent with the experi- 
mentally observed broadening of the NQR line width j| . 
Moreover, we argue that the system undergoes a quan- 
tum phase transition at a critical (3 C , where the ground 
state changes from being unpolarized for [3 < [3 C to a po- 



larized one for (3 > (3 C . Thus, the chain SC state is in 
many aspects similar to that of an s-wave superconduc- 
tor. Finally, we show that in the presence of a second 
impurity, interference effects due to multiple scattering 
processes lead to a splitting of the resonances. 
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ported by DOE at Los Alamos. 
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